arXiv:1501.04844vl [astro-ph.GA] 20 Jan 2015 


Draft version January 21, 2015 

Preprint typeset using style emulateapj v. 12/16/11 


DUST HEATING BY LOW-MASS STARS IN MASSIVE GALAXIES AT Z < 1 

M. Kajisawa^’^, T. Morishita®’®, Y. Taniguchi^, M. A. R. Kobayashi^, T. Ichikawa^, Y. Fukui'^ 

Draft version January 21, 2015 

ABSTRAGT 

Using the Hubble Space Telescope/'Wide Eield Gamera 3 imaging data and multi-wavelength photo¬ 
metric catalog, we investigated the dust temperature of passively evolving and star-forming galaxies at 
0.2 < z < 1.0 in the CANDELS helds. We estimated the stellar radiation field by low-mass stars from 
the stellar mass and surface brightness profile of these galaxies and then calculated their steady-state 
dust temperature. At first, we tested our method using nearby early-type galaxies with the deep EIR 
data by the Herschel Virgo cluster survey and confirmed that the estimated dust temperatures are 
consistent with the observed temperatures within the uncertainty. We then applied the method to 
galaxies at 0.2 < z < 1.0, and found that most of passively evolving galaxies with Mgtar > IO^^Mq 
have a relatively high dust temperature of Tjust > 20 K, for which the formation efficiency of molecular 
hydrogen on the surface of dust grains in the diffuse ISM is expected to be very low from the labora¬ 
tory experiments. The fraction of passively evolving galaxies strongly depends on the expected dust 
temperature at all redshifts and increases rapidly with increasing the temperature around T^nst 20 
K. These results suggest that the dust heating by low-mass stars in massive galaxies plays an impor¬ 
tant role for the continuation of their passive evolution, because the lack of the shielding effect of the 
molecular hydrogen on the UV radiation can prevent the gas cooling and formation of new stars. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: star formation 


1. INTRODUCTION 

Galaxies are generally divided into two populations, 
namely, star-forming galaxies and passivel y evolving 
galax ies with little current star formation ('e.g.. lBell et al.l 
1200411 . At redshift z < 1, most of massive galaxies 
are passively evolving, while lo w-mass galaxies tend to 
actively form new stars (e.g., iKauffmann et all I2003I 
iBundv et al.l 1200^ . Then massive galaxies are consid¬ 
ered to have finished their star formation in the early 
universe, and the star formation can be seen only in grad¬ 
ually less massive galaxies as time passes, the so-calle d 
down-sizing evolution of galaxies (jGowie et al.l 1199611 . 
Since cold gas is expected to successively accrete onto 
massive galaxies in the structure formation model of the 
AGDM universe, such down-sizing effect suggests that 
some quenching mechanisms for the star formation con¬ 
tinue to work in massive galaxies at z < 1. As a possible 
such mechanism, radio-mode feedback, in which the ki¬ 
netic energy of radio jets from active galactic nuclei heats 
the gas in the dark matter halo and prevents its radiative 
cooli ng, has been extensively discussed in many studies 
(e.g., lFabian|[2012D . 

On the other hand, several observational studies re¬ 
ported that the surface stellar mass density of galax¬ 
ies is better correlated with their star formation activity 
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and color tha n stellar mass itself at both high and low 
redshifts (e.g., Kauffmann et al.ll20Ti^ lEranx et al.ll2008l 
I Williams et al. 12010(1 . In fact, passively evolving mas¬ 
sive elliptical galaxies tend to have higher surface mass 
density than star-forming spir al galaxies in the present 
universe ('e.g.. lShen et al.ll2003D . The surface stellar mass 
density is considered to be related closely with the vol¬ 
ume number density of stars within the galaxy, in par¬ 
ticular that of long-lived low-mass stars which domi¬ 
nate the total stellar mass of the galaxy. In this study, 
we relate the surface stellar mass density to the equi¬ 
librium temperature of dust grains heated by low-mass 
stars in galaxies, and discuss its relationship with their 
star formation activity. The formation of molecular hy¬ 
drogen is c onsidered to occur on the surface of dust 
grains le.g.. lGould fc SalDetei([T9^ . and the dust tem¬ 
perature can affec t its formation in the diffuse interstel¬ 
lar medium (e.g., iKatz et ^11999(1 . If the dust grains 
have been heated up to sufficiently high temperatures in 
massive galaxies, the formation of molecular hydrogen 
would continue to be suppressed for a long time, and 
the lack of the shielding effect of the molecular hydro¬ 
gen on the UV radiation would prevent the gas cool¬ 
ing and formation of new stars. In order to examine 
whether the stellar radiation field by low-mass stars is 
sufficiently strong or not, we calculate the stellar radia¬ 
tion field from stellar mass and surface brightness profile 
and then estimate the steady-state dust temperature for 
galaxies at 0.2 < z < 1.0 in the CANDELS survey fields 
(|Grogin et al.l[2Mlll . 

Section 2 describes the sample selection and the 
method to estimate the stellar radiation field and the 
dust temperature. In Section 3, we check our method by 
comparing with the direct measurements of the dust tem¬ 
perature for nearby early-type galaxies by the Herschel 
satellite. We present the results for our main sample of 
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fields. The dotted-dashed line shows the color criteria for the pas¬ 
sively evolving galaxies. Red dots show galaxies with a Sersic index 
of n > 2.5 and blue ones represent galaxies with n < 2.5. 


galaxies at 0.2 < z < 1.0 and discuss them in Section 
4. Throughout this paper, magnitudes are given in the 
AB system. We adopt a flat universe with llmatter = 0.3, 
Oa = 0.7, and Hq = 70 km s“^ Mpc“^. 

2. SAMPLE AND ANALYSIS 

In this study, we used the 3D-HST WFC3-selected pho¬ 
tometric catalog version 4.10 and the WFC3 imaging 
data for the five CANDELS/3D-HST fields (GOODS- 
North and South, AEGIS, COS MOS, and UKIDSS UPS) 
released by the 3D-HST team (|Skelton et al.l[20T^ . The 
source detection was performed on the combined WFC3 
images made from the Jfi 25 Wj and i/pieoW" 

bands data. The 3D-HST team performed the SED fit¬ 
ting and estimated the photometric redshift, rest-frame 
colors, and stellar mass using the photometric data from 
{7-band to Spitzer/IJiAC 8.0 pm band. Chabrier Ini¬ 
tial Mass Function (IMF) (|Chabrieii 1200311 was assumed 
in their calculation of stellar mass of galaxies. At first, 
we selected extended sources with Tlpieow < 25 mag 
and “useJlag” = 1 , which ensures a quality of the pho - 
tometry and SED fitting procedure (|Skelton et al.ll201^ . 
from the catalog to achieve a high completeness for the 
surface brightness profile fitting with the WFC3 images 
described below. We then made a sample of galaxies 
with Mstar > 10 ®Mq at 0.2 < ZpPot < 1-0. The magni¬ 
tude limit of i7pi6ow < 25 is sufficient to sample galaxies 
with Mstar > 1 O®M 0 up to z ^ 1 . 

We divided the sample galaxies into the passively 
evolving and star-forming populations on the rest-frame 
U — V vs. V — J two-colo r diagr am as seen in Figure 
[T1 following IWilliams et al.l (1200911 . The color selection 
criteria for passively evolving galaxies we used are as fol¬ 


lows, 

U-V >0.88x{V-J)+0.59kV-J < 1.6kU-V> 1.4. 

( 1 ) 

The other galaxies are classified as star-forming galaxies. 
We a pplied this selection f or galaxies at 0.2 < z < 1.0, 
while I Williams et al.l (|2009ll changed the criteria slightly 
depending on redshift. We also used the slightly more 
stringent criterion of the U — V color to prevent the 
contamination from star-forming galaxies. We confirmed 
that the small differences in the color criteria do not af¬ 
fect our results at all. 

In order to estimate the dust temperature from the 
stellar radiation field by low-mass stars, we analyzed the 
surface brightness of our sample galaxies by ourselves, us¬ 
ing the WFC3 i7pi6ow-band images. The i7pi60W band 
samples the rest-frame 8000-13000 A light for galax¬ 
ies at 0.2 < z < 1.0. We are interested in the radia¬ 
tion field made by low-mass stars, whose light dominates 
the near-infrared radiation. Therefore the i7pi60W"band 
data are suitable for our purpose. The morphological 
K-correction within these rest-frame wavelengths is ex¬ 
pected to be very small and it does not affect the estimate 
of the dust temperature as s hown in the next sec tion. 
Following our previous study (jMorishita et al.ll2014[l. we 
used the publicly available code GALFIT (iPeng et al] 
l 2002 f) to fit the su rface brightn ess profiles of galaxies with 
the Sersic profile (lSersidll9^ . GALFIT provided us the 
best-fit Sersic index n, half-light semi-major radius Ue, 
and half-light semi-minor radius be- As shown in previ¬ 
ous studies at z < 1, galaxies with n > 2.5 tend to be 
passively evolving, while most galaxies with n < 2.5 are 
star-forming (Figured)). In total, we used 9281 galaxies 
at 0.2 < z < 1.0 with the measured surface brightness 
parameters over a total area of 896.3 arcmin®. Out of 
9281 sample galaxies, 2409 galaxies have n > 2.5, while 
there are 6872 galaxies with n < 2.5. 

We inferred the 3-dimensional distribution of stellar 
mass from the 2 -dimensional surface brightness profile 
mentioned above, in order to estimate the stellar ra¬ 
diation field. For simplicity, we assume a constant 
Mstar/Apieow ratio within a galaxy, which means that 
the profile of the surface stellar mass density has the 
same shape as the Lfpieow-band surface brightness pro¬ 
file. Under this assumption, we can slightly underes¬ 
timate the stellar mass of the bulge and overestimate 
that of the disk for galaxies with a color gradient such 
as disk galaxies with red bulge and blue disk. As a re¬ 
sult, we can underestimate the dust temperature at a 
small radius of the galaxy and overestimate at a large 
radius. However, the effect is relatively small because of 
the relatively small variance of the stellar M/L ratio in 
the rest-frame ^ 10000 A, and in fact our results do not 
change significantly if we measure the dust temperature 
at the different radii such as 0.5 x re and 2.0 x re instead 
of re. From the distribution of the axial ratio bg/ae seen 
in Figured we expect that most galaxies with n > 2.5 
have spheroidal shapes, while galaxies with n < 2.5 tend 
to have thin disk shapes. For the distribution of stars in 
galaxies with n > 2.5, we used a sphe rically symmetric 
Hernquist-like profile (|Hernauistlll990[ ). 


http: //3dhst.research.yale.edu/Data.php 
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Fig. 2.— Distribution of the observed axial ratio be/ue for our 
sample galaxies with the different Sersic indices. 


where and 7 are fitted so that the 2 -dimensional pro¬ 
jection of this profile is matched to the circularized Sersic 
profile with the observed Sersic index n and half-light 
radius ^ is normalization constant de¬ 

termined by requiring that the integration of the profile 
is equal to the total stellar mass of each galaxy. For 
galaxies with n < 2.5, we used the following thin disk 
profile, 


p(r,C) = 2 lexp 


677 


l/r, 


exp 


-bi 


O.lm 


( 3 ) 

The profile along the disk radius is the observed Sersic 
profile and that along the disk height is assumed to be 
the exponential profile. Note that we adopted the semi¬ 
major axis Ue as a radius of the disk to take account of 
the inclination effect. We also assumed a scale height of 
0.1 X Oe, which is a typical value for local disk galaxies. 
Even if we adopt O.OSoe or 0.2ae as a scale height instead 
of O.lttp, the estimated dust temperature changes only by 
- 2-3 % (i.e., 0.3-1.0 K). 

We then calculated the stellar radiation field as 


U{r) = 


M., 


Lho\ 


P{ro) 


1 


47 r|r — rope 


and 


C/(r, 0 ) = 




Tbol 


p(fo,Co) 


1 


47 r|r — rope 


dV (4) 


dV (5) 


for galaxies with n > 2.5 and n < 2.5, respectively. For 
the thin disk model, the radiation field at various radii 
on the disk plane (i.e., at C = 0) was calculated. We 
here used the bolometric stellar mass-to-luminosity ratio 
Tfstar/Tboi of a single 3 Gyr burst model with an age of 
6 Gyr from the GALAXE V population synthesis library 
((Bruzual &: Charl^ 1200311 in order to consider the radi¬ 
ation from long-lived, low-mass stars. As shown in the 


next section, a change of the model age or duration of 
the star formation does not significantly affect our results 
as long as more than 1 Gyr passed since the star forma¬ 
tion stopped. While our sample galaxies are expected to 
have various star formation histories and stellar ages, we 
simply use the single 3 Gyr burst model with an age of 6 
Gyr, because we aim to investigate how the dust temper¬ 
ature is maintained only by low-mass stars after the star 
formation stopped rather than estimate the exact tem¬ 
perature of galaxies at the observed epoch. Therefore, 
for star-forming galaxies, we infer the dust temperature 
in the case that the star formation was quenched by some 
mechanism such as supernova feedback immediately af¬ 
ter the observed epoch and then short-lived massive stars 
died. In other words, we aim to confirm whether the 
dust temperature continues to be high enough to sup¬ 
press the production of the molecular hydrogen and to 
maintain the passive evolution if the star formation has 
once stoppe d in the star- f ormin g galaxies. 

Following iGroves et al.l (I2012I1 . we simply adopted the 
following assumptions to estimate the steady-state dust 
temperature. (1) The dust emission is expressed by the 
modified black body. (2) The stellar radiation field is not 
significantly affected by the dust absorption (i.e., opti¬ 
cally thin). And, (3) the cooling by the dust emission 
balances the heating by the stellar radiation. Although 
star-forming galaxies are expected to have some amount 
of dust extinction, we here consider the contribution only 
from low-mass stars, which radiate their energy mainly 
at the NIR wavelength, and estimate the dust tempera¬ 
ture expected after massive stars died for these galaxies 
as mentioned above. Therefore the assumption (2) is rea¬ 
sonable for our purpose. Under these as sumptions, th e 
dust temperature is expressed as follows (|Drainell2011h . 


Tdust = 



-TT^C 

(Qabs) * 


[60r(4-b/3)C(4+/3)<T 

Qq 


i/(4-e/3) 




where h, feu, and a are the Planck constant, Boltz¬ 
mann constant, and Stephen-Boltzmann constant, re¬ 
spectively. Qo is the dust absorption cross-section at a 
reference wavelength of Aq = 100 /im and (Qabs)» is the 
spectrum-averaged absorption cross-section. The emis- 
siv ity slope of = 2 and the Milky Way dust model 
by lWeingartner fc Drainel (j2001l l with Ry = 3.1 were as¬ 
sumed. We calculated (Qabs)* /Qo 209 from the Milky 
Way dust model and the stellar SED of the single 3 Gyr 
burst model with an age of 6 Gyr. Since the stellar radi¬ 
ation field becomes stronger at smaller radii in a galaxy, 
the estimated dust temperature increases with decreas¬ 
ing radius. In order to ensure the dust temperature is 
higher than a giyen yalue oyer a significant volume of the 
galaxy, we adopt the temperature at a radius of r = rg. 
Even if we use the temperature at a different radius such 
as 0.5 X Te or 2.0 x re, the distribution of the dust tem¬ 
perature only slightly shifts to higher or lower values as 
shown in the next section and the results in this paper 
do not change. 

We also checked the effects of our simplification for 
the 3-dimensional distribution of stellar mass by com¬ 
paring the dust temperatures estimated with the spher¬ 
ically symmetric model and the thin disk model of the 
same objects. In Figure [3l we compare the dust temper- 
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Expected dust temperature (thin disk model) 

Fig. 3.— Comparison between the expected dust temperatures 
estimated with the thin disk model and the spherically symmetric 
model of the 3-dimensional distribution of stars for galaxies with 
n = 2.0-3.0. 

atures calculated with the spherical model and the thin 
disk model for galaxies with n = 2.0-3.0. Our spheri¬ 
cal model predicts slightly higher temperatures than the 
thin disk model, and the mean and dispersion of the dif¬ 
ferences between these models are AT = 1.9±1.2 K. The 
difference depends on the axial ratio, and the difference is 
very small (AT < 1 K) for objects with foe/oe ~ 1, while 
it becomes larger (AT ~ 4 K) for those with foe/oe ~ 
0.2-0.4. Since the differences between the two extreme 
cases (i.e., spherically symmetry and thin disk) is rel¬ 
atively small, we expect that our simplification for the 
distribution of stars does not significantly affect our re¬ 
sults. In fact, we obtain nearly the same results if we 
divide galaxies into the spheroidal shapes and thin disks 
at n = 2.0 or 3.0 instead of n = 2.5. 

3. COMPARISON WITH DIRECT MEASUREMENTS 

In this section, we check the method described in 
the previous section, using the direct measurements 
of the dust temperature of l ocal early-type galaxies. 
We us ed the measurements bv Idi Serego Alighieri et ahl 
(IMl for early-type galaxie s in the Virgo cluster. 
Idi Serego Alighieri et al.1 (I2013D fitted the Herschel 100, 
160, 250, 350, and 500 /rm-ba nds photometry fro m the 
Herschel Virgo cluster survey (jPavies et al.ll2?)T^ with 
the modified black body with /3 = 2 to measure the tem¬ 
perature of cold dust. Their deep 5-bands data from 100 
to 500 /im enable the estimation of the dust temperature 
with relatively high accuracy. Since their sample has a 
wide range of stellar mass, we can expect a wide range 
of stellar surface density and dust temperature, which is 
essential for our purpos e, although the sample is li mited 
to the cluster galaxies. Idi Serego Alighieri et al.l (l2013f ) 
also estimated the stellar mass of the sample galaxies 
from their optical-NIR colors. The estimate of the stellar 
M /L ratio was based on the population synthesis model 
by iBruzual fc Charl^ (|2003D and the Chabrier IMF was 



Observed dust temperature (K) 

Fi g. 4.— Comparison between th e observed dust temperature 
from Idi Serego Alighieri et aLl <120131) and the expected tempera¬ 
ture which we estimated from the stellar mass and surface bright¬ 
ness profile for early-type galaxies in the Virgo cluster. The upper 
(bottom) panel shows the expected temperature for which we used 
the SED model of the single 3 Gyr burst model with an age of 6 Gyr 
(10 Gyr) in the calculation. Gircles show galaxies with the Sersic 
index of n > 2.5, while triangles represent those with n < 2.5. The 
errors in the expecte d temperatures are based on th e uncertainty 
in stellar mass from Idi Serego Alighieri et ahl 1)20131) and that in 
the radius for which we assume 30 % error taking account of the 
systematic uncertainty due to the errors in the background esti¬ 
mate. 


assumed. We us ed the dust temperature and ste llar mass 
from Table 4 of idi Serego Alighieri et al] (|2013D . 

In order to measure the surface brightness profile of 
these galaxies, we used z- band images froin the Sloan 
Digit al Sky Survey (SDSS; lYork et S] 120001: lAhn et ahl 
I2014D . These data have a pixel scale of 0.4 arcsec/pixel 
and the PSF FWHM of ^ 1.5 arcsec. We performed 
the same procedure of the surface brightness fitting with 
GALFIT to determine the best-fit Sersic index n, half- 
light semi-major radius Ue, and half-light semi-minor ra- 
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dius be- Using the fitted parameters of the Sersic profile 
from the SPSS z-band data an d the stellar mass from 
tdi Serego Alighieri et alj (I2013D . we then estimated the 
dust temperature at a radius of r = Tg as described in 
the previous section. We excluded the objects for which 
the dust temp erature cannot be deterrnined b y the FIR 
SED fitting of fdi Serego Alighieri et all (l2013f l or the ef¬ 
fective radius cannot be reliably derived in the surface 
brightness fitting, for example, due to their faintness on 
the SDSS images. 

The upper panel of Figure |4] shows a compar- 
ison between the obse r ved dust temperature from 
Idi Serego Alighieri et al.l (I2013D and the expected dust 
temperature we derived from the stellar mass and surface 
brightness profile. It is seen that the expected dust tem¬ 
perature correlates relatively well with the temperature 
estimated from the FIR observations, although there is 
some dispersion. The mean and dispersion of the differ¬ 
ences between the observed and expected temperatures 
are AT = —0.4 ±3.3 K. We also estimated the dust tem¬ 
perature using the single 3 Gyr burst model with an age 
of 10 Gyr as the star formation history instead of that 
with an age of 6 Gyr, because the 10 Gyr age can be 
more suitable for these early-type galaxies in the clus¬ 
ter at z ^ 0. The result is shown in the bottom panel 
of Figure m The expected temperatures become slightly 
lower values than those in the case with the 6 Gyr age. 
In this case, the difference between the observed and ex¬ 
pected temperatures is AT = —1.9 ± 3.1 K. The mean 
difference between the temperatures estimated from the 
model SEDs with the 6 Gyr and 10 Gyr ages is ^ 1.5 K, 
which is smaller than the dispersion around the mean. 
As mentioned in the previous section, the difference in 
the assumed model age does not affect so strongly the 
estimated dust temperature. 

We also checked the expected dust temperatures at the 
different radii, namely r = 0.5 x rg and r = 2.0 x Cg in the 
galaxies. The differences between the observed tempera¬ 
ture and the expected ones with the 6 Gyr age model are 
AT = 2.8 ± 4.1 K for r = 0.5 x Cg and AT = -3.9 ± 2.8 
K for r = 2.0 x rg. The expected dust temperature be¬ 
comes higher or lower by ~ 3-4 K on average, when we 
use the temperature at r = 0.5 x rg or r = 2.0 x rg. 
The results in the next section do not change even if 
we use the different radii. The expected dust tempera¬ 
ture at r = rg seems to show a better agreement with 
the observed temperature than those at the smaller and 
larger radii, which may be consistent with the relatively 
compact morphology of the FIR dust e mission in these 
galaxies (jdi Serego Alighieri et al.|[2013ll . 

Furthermore, we performed the same analysis with J- 
band images from the 2-Micron All Sky Survey (2MASS; 
iSkrutskie et al.l [20061 1 in order to examine the effect of 
the morphological K-correction on the estimate of the 
dust temperature. The pixel scale of these data is 1.0 
arcsec/pixel and the PSF FWHM is ^ 3 arcsec. The 
difference between the observed and expected tempera¬ 
tures is AT = 1.5 ± 4.1 K. The difference between those 
with the z-band and J-band images is ^ 2 K on average, 
which is considered to be relatively small if we take ac¬ 
count of the dispersion around the mean and the different 
depth and spatial resolution between the SDSS z-band 
and 2MASS J-band data. 

In summary, the dust temperature estimated from the 


stellar mass and surface brightness profile agrees rela¬ 
tively well with the observed temperature for early-type 
galaxies in the Virgo cluster with a dispersion of ~ 3- 
4 K. There can be also the systematics of ~ 3-4 K in 
the estimated temperature, depending on the assumed 
star formation history in the SED model, the wavelength 
where the surface brightness is measured, and the radius 
used in the calculation of the dust temperature. Keep¬ 
ing in mind these uncertainties, we discuss the expected 
dust temperature for passively evolving and star-forming 
galaxies at 0.2 < z < 1.0 in the following section. 

4. RESULT AND DISCUSSION 

We show the distribution of the surface stellar mass 
density and expected dust temperature as a function 
of stellar mass for passively evolving and star-forming 
galaxies at 0. 2 < z < 1.0 in Figure El As previously 
reported fe.g. iKauffmann et al.ll200^ iFranx et al.l 1200^ 
IWilliams et al.ll2010tl . most galaxies with a high surface 
stellar mass density of Egtar ^ 10®MQ/kpc^ are passively 
evolving ones especially at Mgtar > IO^^Mq. Similarly, 
most of passively evolving galaxies show a relatively high 
dust temperature of Tdust > 20 K at Mgtar > 10 ^°Mq. 
The expected dust temperature at smaller radii is higher 
because the stellar radiation increases with decreasing 
radius. Since we adopted the temperature at a radius of 
r = Te in the figure, the dust temperature is expected to 
be higher than 20 K over a significant volume in most of 
massive passively evolving galaxies. 

Since the energy density by the radiation from a star 
with a luminosity T at a distance of r is expressed by 
T/47rr^c, the stellar radiation field by all stars at a given 
point in the galaxy is proportional to Mstar/R^, where R 
is the size of the galaxy. Therefore, the stellar radiation 
field in galaxies is directly related with their surface stel¬ 
lar mass density, and the close relationship between the 
dust temperature and the surface stellar mass density is 
naturally expected. 

Figure [6] shows the rest-frame U — V color as a func¬ 
tion of the expected dust temperature for galaxies with 
Mstar > 10l°Mo in the different redshift bins. The rest- 
frame color becomes redder with increasing the dust tem¬ 
perature at all redshifts, and passively evolving galax¬ 
ies on the red sequence dominate at Tdust ^ 20-25 K. 
The transition temperature between star-forming and 
passively evolving populations does not seem to signif¬ 
icantly evolve with time. We show the fraction of pas¬ 
sively evolving galaxies as a function of the dust temper¬ 
ature in Figure [71 The fraction strongly depends on the 
dust temperature and increases rapidly with increasing 
the temperature around Tdust 20 K. 

Recent observations of nearby early-type galaxies 
by the Herschel satellite revealed that massive early- 
type galaxies in fact tend to have a rel atively high 
dust te mperature o f Td„st > 20 K (e.g., [Sm ith et al. 
20121 : lAuld et al.l 2013 : Idi Serego Aliriiieri et al. 

2013t lAmblard et al.l 2014 1. [Smith et al.r()2012[l and 
Auld et al.l (j201.'^ reported that early-type galaxies 
show systematically higher dust temperature s than 
late-type galaxies. Idi Serego Alighieri et al.l (j2013f l 
suggested that massive early-type galaxies tend to have 
higher dust temper ature than low-mass early-type ones. 
iGroves et al.l (1201 2f) also reported that the bulge of M31, 
which has no young stellar population, shows a high 
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Fig. 5. — Distribution of surface stellar mass density (upper 
panel) and the expected dust temperature at a radius of r = 2.0 X re 
(lower panel) of galaxies at 0.2 < z < 1.0 in the five CANDELS 
fields as a function of stellar mass. Color of the symbols shows 
passively evolving (red) and star-forming (light blue) populations. 
Circles show galaxies with the Sersic index of n > 2.5, while trian¬ 
gles represent those with n < 2.5. The surface stellar mass density 
is calculated by Afstar/27rrg, which corresponds to the averaged 
value within a half-light radius. 


dust temperature of 20-35 K, while the temperature in 
the star-forming disk is ^ 17 K. They found that the 
temperature profile of the M31 bulge can be explained 
by the stellar radiation field from its old stars. These 
studies are considered to be consistent with our results. 

The dust temperature Tdust ^ 20 K corresponds to 
a temperature above which the formation efficiency of 
the molecular hydrogen in the diffuse ISM is expected 
to be very low from t he laboratory exp er iments (e.g., 
Pirronello et al.l 119991: iKatz et al.l llOoiH iPerets et alJ 
2005t iLe Petit et al.ll2009D . The molecular hydrogen can¬ 
not be efficiently formed in the gas phase, and its forma¬ 
tion is considered to occ ur on the surface of dust grains, 
which work as catalysts (|Gould fc SalDeteilll963|) . How- 



Fig. 6. — The rest-frame U — V color of galaxies with Mstar > 
as a function of the expected dust temperature for the 
different redshift bins. Symbols are the same as Figure [S] 
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Fig. 7.— Fraction of passively evolving galaxies in galaxies with 
Mstsur > (solid circle) and Mstar > 10 ®Mq (open square) 

at 0.2 < 2 < 1.0 as a function of the expected dust temperature. 
Data points for the temperature bins which include more than 
10 passively evolving galaxies are plotted. The total number of 
galaxies with Mgtar > 10^® Mq (shaded histogram) and M^tar > 
10 ®Mq (open histogram) is also shown in the upper panel. 

ever, the hydrogen atom that sticks to the surface of 
the grain quickly desorbs and cannot remain on the sur¬ 
face to encounter another atom and form a molecule at 
Tdust ^ 20 K, although this may not be the case for more 
dens e photo-dissociation reg i ons of star-fo rming galaxies 
(e.g., iLe Bourlot et al.ll20ia lTielensl[201,'lD . Thus if low- 
mass stars in massive galaxies has heated up the dust 
grains to Tdust > 20 K, the formation of the molecular 
hydrogen could have been suppressed in the diffuse ISM 
of these galaxies. 

Such a suppression of the formation of the molecular 
hydrogen is expected to play an important role for the 
passive evolution of massive galaxies. If we consider that 
the gas is accreting and cooling from higher tempera¬ 
tures after the star formation has once been stopped by, 
for example, intense supernova feedback, the HI gas is 
expected to settle to the equilibrium temperature (e.g., 
Tgas 70 K for the cold ne utral medium in the Milky 
Way; iHeiles fc Trolandll2003D determined from a balance 
between the heating mainly by photoelectron from dust 
grains and cosmic-ray ionization and the radiative cool¬ 
ing by the emission lines such as [Cll] 158 /rm and [Ol] 
63 /im (lDrainel[2011l ). If the formation of the molecular 
hydrogen proceeds in the cold HI gas, the self-shielding 
effect of the molecular hydrogen on the UV radiation 
would suppress the heating, which leads to further cool¬ 
ing of the gas and the star formation. On the other 
hand, if the formation of the molecular hydrogen in the 


diffuse ISM is suppressed due to the high dust tempera¬ 
ture, the HI gas would keep the equilibrium temperature. 
Therefore the dust heating by long-lived, low-mass stars 
in massive galaxies can cause the continuation of their 
passive evolution. Recent observational studies of the HI 
gas in local early-type galaxies found that a significant 
fraction of massive early-type galaxie s have a moderat e 
amount of HI gas (e.g., IO^-IO^^Mq; iSerra et aI1l20I2ll . 
The low star-formation activity of such galaxies may be 
related with their high dust temperature. 

On the other hand, several studies with radio obser¬ 
vations also detected the CO lines in a significant frac¬ 
tion of nearby earl y-type galaxies ('e.g.. lWelch et al.ll201Ct 
lYoung et al.ll20Ill) . The morphology of the detected CO 
lines tends to be centrally concentrated and those early- 
type galaxies with the CO detection s how the star for¬ 
mation activity at their central region (jSerra et al.ll2ni2l : 
lYoung et al.l 120141 ) . In the case that the cold gas can 
be concentrated into the center of the galaxy, for exam¬ 
ple, by the galaxy interactions, the molecular hydrogen 
may be formed and star formation may occur irrespec¬ 
tive of the dust temperature, because the high density 
of the gas and dust can lead to a enough shielding for 
the UV radiation by the dust absorption and the gas can 
be cooled to form the molecular hydrogen. Otherwise, 
some of the detected molecular gas may be explained by 
the externa l origin, namely, t h e cold gas ac c retion into 
the galaxy (jYoung et al.ll2oi^ . lYoung et al.l (j20lH) also 
found that the molecular gas is detected preferentially 
in the fast rotator early-type galaxies, while there is few 
slow rotators with the CO detection. Therefore the mass 
assembly history may be related with the existence of the 
molecular gas in early-type galaxies. Since the surface 
stellar mass density of the fast r otators tends to be lowe r 
than that of the slow rotators (jCapDellari et al.ll20I3D . 
the dust temperature may cause the difference between 
these two populations in the CO detection rate. 

It is noted that low-mass passively evolving galaxies 
have the similar surface mass density with star-forming 
galaxies, and therefore their dust temperature is ex¬ 
pected to be correspondingly low (see Figure [5|). The 
passive evolution of these low-mass galaxies may be ex¬ 
plained by other mecha nisms such as env ironmental or 
satellite quenching (e.g.. lPeng et al.l[2012l) . 

While we assumed the SED of the old stellar pop¬ 
ulation to investigate the effect of the radiation from 
low-mass stars, there are young massive stars in star¬ 
forming galaxies. Although these massive stars are much 
luminous than low-mass stars and heat the surround¬ 
ing dust, the ultraviolet light from these stars is effi¬ 
ciently absorbed by the surrounding dust and molecu¬ 
lar hydrogen (self-shielding). Therefore the temperature 
of the cold dust in disks of spiral galaxies tends to be 
relativelY low (Th, 17 K; iTabatabaei fc BerkhuiisenI 
120101: iClemens et al.l[^013[) . while there is also the warm 
dust component. On the other hand, such self-shielding 
effect for the near-infrared radiation from low-mass stars 
is expected to be much lower. Thus the radiation from 
low-mass stars may play a role to heat the dust globally 
in star-forming galaxies when low-mass stars are suffi¬ 
ciently accumulated as the star formation proceeds. 

The dust heating by low-mass stars can contribute to 
the passive evolu tion of massive compact galaxies found 
at z ^ 2 (e.g., iTruiillo et al.l 120(371 : Ivan der Wei et al] 
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[Ml . Since a very strong stellar radiation field is ex¬ 
pected within such compact massive galaxies, the dust 
temperature is probably high in these objects. The effect 
of the high dust temperature may also be consistent with 
the i nside-out scenario of the massive galaxy for mation 
le.g.. Ivan Dokknm et al.ll201fll : iPatel et al.l[2013fl . The 
star formation outside a compact core could be permit¬ 
ted if the radiation field is relatively low at outer regions, 
while the formation of molecular hydrogen is strongly 
suppressed within the core. The direct measurements of 
the dust temperature and its profile for massive galaxies 
at high redshift by ALMA will be important to verify 
these scenarios. 

We would like to thank the referee for many invalu¬ 
able suggestions. This work is based on observations 
taken by the 3D-HST Treasury Program (GO 12177 
and 12328) with the NASA/ESA HST, which is op¬ 
erated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS5-26555. 
Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, and the U.S. Department 
of Energy Office of Science. The SDSS-III web site is 
http://www.sdss3.org/, SDSS-III is managed by the 
Astrophysicai Research Consortium for the Participat¬ 
ing Institutions of the SDSS-III Collaboration including 


Ahn, C. P., AlexandrofI, R., Allende Prieto, C., et al. 2014, ApJS, 
211, 17 

Amblard, A., Riguccini, L., Temi, P., et al. 2014, ApJ, 783, 135 
Auld, R., Bianchi, S., Smith, M. W. L., et al. 2013, MNRAS, 428, 
1880 

Bell, E. F., Wolf, C., Meisenheimer, K., et al. 2004, ApJ, 608, 752 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bundy, K., Ellis, R. S., Conselice, C. J., et al. 2006, ApJ, 651, 120 
Cappellari, M., McDermid, R. M., Alatalo, K., et al. 2013, 
MNRAS, 432, 1862 
Chabrier, G. 2003, PASP, 115, 763 

Clemens, M. S., Negrello, M., De Zotti, G., et al. 2013, MNRAS, 
433, 695 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 
112, 839 

Davies, J. L, Bianchi, S., Cortese, L., et al. 2012, MNRAS, 419, 
3505 

di Serego Alighieri, S., Bianchi, S., Pappalardo, C., et al. 2013, 
A&A, 552, AA8 

Draine, B. T. 2011, Physics of the Interstellar and Intergalactic 
Medium by Bruce T. Draine. Princeton University Press, 2011. 
Fabian, A. C. 2012, ARA&A, 50, 455 

Franx, M., van Dokkum, P. G., Schreiber, N. M. F., et al. 2008, 
ApJ, 688, 770 

Gould, R. J., &; Salpeter, E. E. 1963, ApJ, 138, 393 
Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 
197, 35 

Groves, B., Krause, O., Sandstrom, K., et al. 2012, MNRAS, 426, 
892 

Heiles, C., & Troland, T. H. 2003, ApJ, 586, 1067 
Hernquist, L. 1990, ApJ, 356, 359 

Katz, N., Furman, L, Blham, O., Pirronello, V., & Vidali, G. 

1999, ApJ, 522, 305 

Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, 
MNRAS, 341, 54 

Kauffmann, G., Heckman, T. M., De Lucia, G., et al. 2006, 
MNRAS, 367, 1394 

Le Bourlot, J., Le Petit, F., Pinto, C., Roueff, E., &; Roy, F. 2012, 
A&A, 541, AA76 

Le Petit, F., Barzel, B., Biham, O., Roueff, E., & Le Bourlot, J. 
2009, A&A, 505, 1153 

Morishita, T., Ichikawa, T., & Kajisawa, M. 2014, ApJ, 785, 18 


the University of Arizona, the Brazilian Participation 
Group, Brookhaven National Laboratory, Carnegie Mel¬ 
lon University, University of Florida, the French Partici¬ 
pation Group, the German Participation Group, Harvard 
University, the Institute de Astrofisica de Canarias, the 
Michigan State/Notre Dame/JINA Participation Group, 
Johns Hopkins University, Lawrence Berkeley National 
Laboratory, Max Planck Institute for Astrophysics, Max 
Planck Institute for Extraterrestrial Physics, New Mex¬ 
ico State University, New York University, Ohio State 
University, Pennsylvania State University, University of 
Portsmouth, Princeton University, the Spanish Partici¬ 
pation Group, University of Tokyo, University of Utah, 
Vanderbilt University, University of Virginia, University 
of Washington, and Yale University. This publication 
makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analy¬ 
sis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. Data analysis were in 
part carried out on common use data analysis computer 
system at the Astronomy Data Center, ADC, of the Na¬ 
tional Astronomical Observatory of Japan. YT acknowl¬ 
edges the financial support from the Japan Society for 
the Promotion of Science (No. 23244031). 


Patel, S. G., van Dokkum, P. G., Franx, M., et al. 2013, ApJ, 

766, 15 

Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 

124, 266 

Peng, Y.-j., Lilly, S. J., Renzini, A., & Carollo, M. 2012, ApJ, 

757, 4 

Perets, H. B., Biham, O., Manico, G., et al. 2005, ApJ, 627, 850 

Pirronello, V., Liu, C., Roser, J. E., & Vidali, G. 1999, A&A, 344, 
681 

Serra, P., Oosterloo, T., Morganti, R., et al. 2012, MNRAS, 422, 
1835 

Sersic, J. L. 1963, Boletin de la Asociacion Argentina de 
Astronomia La Plata Argentina, 6, 41 

Shen, S., Mo, H. J., White, S. D. M., et al. 2003, MNRAS, 343, 
978 

Skelton, R. E., Whitaker, K. E., Momcheva, L G., et al. 2014, 
ApJS, 214, 24 

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 
1163 

Smith, M. W. L., Gomez, H. L., Eales, S. A., et al. 2012, ApJ, 
748, 123 

Tabatabaei, F. S., & Berkhuijsen, E. M. 2010, A&A, 517, AA77 

Tielens, A. G. G. M. 2013, Reviews of Modern Physics, 85, 1021 

Trujillo, L, Conselice, C. J., Bundy, K., et al. 2007, MNRAS, 382, 
109 

van der Wei, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 
788, 28 

van Dokkum, P. G., Whitaker, K. E., Brammer, G., et al. 2010, 
ApJ, 709, 1018 

Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 

Welch, G. A., Sage, L. J., & Young, L. M. 2010, ApJ, 725, 100 

Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & 
Labbe, 1. 2009, ApJ, 691, 1879 

Williams, R. J., Quadri, R. F., Franx, M., et al. 2010, ApJ, 713, 
738 

York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 
120, 1579 

Young, L. M., Bureau, M., Davis, T. A., et al. 2011, MNRAS, 

414, 940 

Young, L. M., Scott, N., Serra, P., et al. 2014, MNRAS, 444, 3408 








